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Abstract 

Let P := {Xi}i>i be a stationary Poisson point process in {Cj}j>i 
be a sequence of i.i.d. random sets in K , and {Yl; t > 0}j>i be i.i.d. 
{0, l}-valued continuous time stationary Markov chains. We define the 
Markov-Poisson-Boolean model Ct := {Y^(X{ + Ci),i > 1}. Ct represents 
the coverage process at time t. We first obtain limit laws for /c-coverage 
of an area at an arbitrary instant. We then obtain the limit laws for the 
£;-coverage seen by a particle as it moves along a one-dimensional path. 
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1 Introduction 



1.1 Motivation 

This paper is motivated by the need to characterize the ability of a sensor network with randomly 
deployed unreliable sensor nodes to track the trajectory of a moving target in the sensor field. Since 
the nodes are randomly deployed, a suitable point process in a ci-dimensional metric space, typically 
$l d or Z d , can be used to describe the location of the sensor nodes. In this paper we will assume 
a homogenous Poisson point process for the distribution of sensor locations. A sensor node can 
detect events or perform measurements over a 'sensing area' or a 'footprint.' The coverage area of 
each sensor is described by a suitable sequence of random (f-dimensional sets. Thus analyzing the 
coverage by the sensor network involves the analysis of an equivalent stochastic coverage process. 
Coverage analysis usually takes the form of obtaining statistics for the fraction of the volume of a 
(i-dimensional set that is covered by one or more sensors, e.g., as in pQ. 

Tracking a moving target by a sensor network involves trajectory estimation from a sequence of 
position estimates of the target. The quality of the trajectory estimate will depend on the parts of 
the trajectory that are 'covered' by k or more sensors with the value of k being determined by the 
estimator used. 

The sensor nodes are assumed to be unreliable in the sense that they toggle between two states — 
'available' and 'not available.' Hence some of the sensors that are covering the target as it is moving 
along the path could become unavailable during the coverage of the target. There are many reasons 
for this. A sensor could be turned off for energy saving or even energy restoration in the battery. 
Alternatively, a sensor may have temporarily failed. It could also be that the radio environment is 
such as to prevent a sensor from communicating with its neighbors, effectively making it unavailable 
for the sensing process. Since the nodes are unreliable, the coverage of the path during the motion 
of the target is a random process in time that is determined by the switching of the sensor states 
that could have covered the path. Thus we need to explicitly model the temporal behavior of the 
induced-coverage process of the path. We model this "on-off" behavior of the sensors by a two 
state Markov chain. 

The aim of this paper is two-fold. We first extend the asymptotic coverage results in pQ to the case 
of general /^-coverage. Asymptotic properties of the covered fraction are derived as the intensity 
of the Poisson point process becomes large and the coverage area of the individual sensors scaled 
down in such a way that the limiting process covers a non-trivial fraction of the operational area. 
The second objective is to show how the dynamics of the on-off process affects the coverage of a 
linearly moving target. 
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1.2 The Markov-Poisson-Boolean Model 

Let P := {Xi, i = 1,2,...} be a stationary Poisson process in of intensity A > 0. The points 
of P can be thought of as locations of sensors in a random sensor network. Let C be a random 
closed set in independent of P and having an arbitrary distribution. Throughout the paper we 
assume that for some fixed r > 0, C C Bq <t , almost surely, where Bq )T is a closed ball of radius r 
centered at the origin. We will also assume that ||C|| > with probability 1, where || • || denotes 
the Lebesgue measure. Let (3 = E(||C||) . Let Cj be i.i.d. copies of C. As in p], the Cj's will be 
called shapes to distinguish them from the sets (Xi + Cj)'s that denote the areas in $t d that are 
covered by the sensors. The coverage process C = {Xi + d,i > 1} is called the Poisson-Boolean 
model [lj. Now let Y = {Y l , t > 0} be a {0, l}-valued continuous time stationary Markov process 
independent of all other random variables. Let {Y?,t > 0},i > 1 be i.i.d. copies of Y. Y? = 1 can 
be interpreted to mean that sensor i is 'on' at time t and is available for sensing and Yf = means 
that it is 'off' at time t and not available for sensing. Define the Markov-Poisson-Boolean model 
Ct := {Y?(Xi + Ci),i > 1}. Ct represents the coverage process by the available sensors at time t. 
Let R C 9f? d be an arbitrary Borel set. R could correspond to the operational area of the sensor 
network. Now, consider a straight line path C in R of length L units, and an object moving along 
it with a constant velocity of c units per second. The object starts moving at time t = 0. For any 
positive integer k, let Ck(t) be the indicator function for the object being 'covered,' to be made 
more precise later, by k or more sensors at time t. The objective of this paper is to characterize 
the random variables Ck(t) an d Jq Ck(t)dt, where T = L/c. In the rest of the paper, we refer to the 
Poisson-Boolean model and the Markov-Poisson-Boolean model as the PB and the MPB models 
respectively. 

1.3 Previous Work 

Coverage with reliable nodes, where nodes are always in the 'available' state, has been elaborated 
in both applied mathematics and sensor network literature. See [T] for a good comprehensive first 
study and [3j for some recent results for the case k = 1. In [ID] . R is a square of side-length I and 
C is a circle of unit area. It is shown for this PB model, that if A is given by 

A = log£ 2 + (A; + l)loglog^ + Q, (1) 

and if ci — > oo as £ — > oo, then R is almost surely asymptotically /c-covered. 

Asymptotic coverage by unreliable sensor networks at an arbitrary epoch has been studied in 0[2]. 
In the analysis of coverage by an unreliable sensor network at an arbitrary epoch, the stationary 
probability of a sensor node not being available essentially 'thins' the original deployment process 
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and a standard analysis with the thinned process applies. However, for applications like target 
tracking or intruder detection, we need to know the behavior of the coverage process during the 
movement of the target. When the sensor nodes are unreliable, a node that was sensing the object 
may switch from being available to becoming unavailable, or vice versa. This implies that the 
coverage of a given point is not independent either in space or in time. Thus we need to consider 
the dynamics of the transitions from availability to non-availability of the sensors in the spatio- 
temporal analysis of the coverage of the path. The coverage of a line by a two-dimensional PB 
model was investigated in [U [5] . 

1.4 Organization of the Paper and Summary of Results 

In Section [21 we characterize /c-coverage in (/-dimensions. Although our eventual interest is the 
characterization of the coverage of a moving point on a path by the MPB model defined earlier, 
it is instructive to first consider the /c-coverage of R C 3? rf , by the MPB model Ct at an arbitrary 
instant. At an arbitrary instant, let (R) be the volume of an arbitrary ci-dimensional set R C 3f? d 
that is not covered by k or more sensors. We obtain a strong law and central limit theorem for 
V/. (R) . The proof techniques are in general similar to those in Chapter 3 of pQ . In the second part 
of Subsection 12.21 we consider the special case when the coverage areas are discs of fixed radius r, 
and obtain a strong law of large numbers for the critical radius required for complete /c-coverage in 
dimension d = 2. 

In Section [3l we consider the MPB model defined earlier, which is the PB model but now with 
unreliable sensors. We analyze the path coverage for a linearly moving target that is in the sensor 
field for T units of time. Without loss of generality, let this interval be (0, T). Let Vk T De the total 
time in (0, T) that the target is not tracked by k or more sensors. For pedagogic convenience we 
consider k = 1 and obtain a strong law and central limit theorem for V\ t T- The proof techniques of 
Section [2] can be combined with those of Section [3] to extend the results to the case of k > 1. We 
have separated /c-coverage from the Markovian on-off dynamics to maintain clarity of exposition. 
Each of these two components operates independently in the computations and the expressions 
involving /c-coverage are at times lengthy. 
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2 A;- Coverage 



2.1 Preliminaries 

For a point x E let Xm( x ) be the indicator variable that x G + C; for exactly m points in P, 



i.e., 



1 if x E Xi + Cj for exactly m of £ = 1, 2, 
otherwise. 



E(% m (x)) = Pr (x € (Xj + Cj) for exactly m points) 
= Pr (for exactly m points X{ £ (x — Cj)) 

= Pr (for exactly m points Xi S Cj) . (2) 

The last equality follows from the stationarity of P. Recall that Cj C Bq t and (3 = E(||C||) . Given 
that N points of P lie within Bq, t , the probability that exactly m of these points cover the origin 
is given by 

N\ ( \ N ~ m ( (3 



K m) \ ||Bo,r||/ \||-Bo,t| 
Since N is Poisson with mean A||-Bor||> we obtain 



mj \ \\Bq iT \\J \I|A),t||/ / TH 

Let R G be a (i-dimensional set. For k > 0, we define the A;- vacancy within R, V k {R), to be the 
c?-dimensional volume of the part covered by at most k — 1 random sets of C, i.e., 



fc-i 

R 



m=0 



V k (R) = V / Xm{x) dx. (4) 



The indicator variable for the ^-vacancy of a point x will be denoted by Vf~ (x), i.e., Vfc(x) := 
Em^o Xm(^)- — Vk(R) will be called the k-coverage of i?. Since R is fixed throughout the 
paper, we will omit the reference to it in the notation and write Vk(R) as 

Some of the early derivations mimic that in [I] and we give it here for the sake of completeness. 
From (l3l) and Fubini's theorem, 



E(V k ) = £ / E( Xm (x))dx=22\\ R \\ ^T— ( 5 ) 



/r ' — ' ml 

m=0 1 m=0 
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We now derive the variance of If a point x is covered by Xj, then x £ Xj + Cj or Xj £ x — Cj. 
Similarly, if x is not covered by Xj then X,- L £ x — Cf where C c = \ C. We use this to first obtain 
the probability that two points x\ and X2 are covered by exactly m and n sensors respectively 
which is then used to obtain the variance of We make the following observations regarding the 
location of X; t relative to points x\ and X2- 

• If Xi covers x\ and X2, then Xj £ B\(xi,X2) ■= (xi — Cj) n (x2 — Cj). Further, has the 
same distribution as where B\{x\ — x^) = (x± — X2 + C) H C. 

• If Xj covers x\ and not X2, then Xj £ B\{x\,X2) := (^i — Cj) n (x2 — Cf) and H-B2II ^ as ^ ne 
same distribution as H-B2II where B2(x± — X2) = (x% — X2 + C) n C c . 

• Similarly, if Xj does not cover xi but covers £2, then Xj £ B\(x\, X2) '■= {x\ — Cf ) (X2 — Cj) 
and H-B3II is equal to H-B3H in distribution where B^{x\ — X2) = (x± — X2 + C c ) n C. 

We will suppress the argument of Bj and B l - unless required. Observe that the Bj defined above 
are mutually disjoint sets. Further, H-B2II an d \\B^\\ will have the same distribution. 

E(xm(xi)Xn(x2)) = Pr {x\ covered by m & X2 covered by n sensors) 

= Pr (xi G Xj + Cj for m sensors and X2 £ Xj + Cj for n sensors) 
= Pr {Xj £ xi — Cj for m sensors and Xj £ X2 — Cj for n sensors) 

mhn 

= P r (l cover x\ & X2; (m — I) cover only x±; (n — I) cover only X2) . 
1=0 

We can proceed as in the derivation of (|2|) and ([3]) and consider a bounded set A' that contains [x\ + 
Bq t ) and (x2+So r ). The 'left' point is designated x\. Let N be the number of points of the Poisson 
point process P of intensity A lying in A'. Then the probability that (m—l) sensors cover x\ only and 
(n—l) cover X2 only and I sensors cover both is Pr (Mi = I, M2 = m — I, M3 = n — I, M4 = N — m — n + l) 
with 

(Mi, M2, M3, M4) ~ Multinomial(X, a\, 02, 03, 04 = 1 — a\ — 02 — 03), 

where 

Since X is Poisson with mean A||A'||, the unconditional probability that (m — I) sensors cover xi 
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but not x 2 , (n — I) cover x 2 but not xi, and I sensors cover both x\ and X2 will be 

^ / (m-O n (n-/) n (^-m-n+0^ _ ( AE ( II g l II ))' _-Ag||Bi I 

n( m _/)!( n _/)!(Ar-m-n + /)! aia2 ° 3 ° 4 J" Z! 

(m-Z)! (n-Z)! 
Hence from above calculations, 

E(x„(x l) X„(x 2 )) - "f ™^ 

z=o ' ^ m '' 

(n-Z)! 



,-2^11*11) V (AEdl^H))' (AEC^r^ 



Z=0 



The last equality is obtained from the identities E(||5i|| + ||B 2 ||) = E(||5i|| + ||B 3 ||) = (3. We then 
have 



fe-i 



Cov(y fc (xi),Vfc(x 2 )) = X] CoV (Xm(Sl),Xn(32)) 



ra,n=0 
fc-1 



= ( E (Xm(a;i)Xn(a;2)) - E(xm(xi)) E(xn(x 2 ))) 



m,n=0 

g -2A/3 ^ ^ c AEri|Bifxi-x 2 )||) ^ (AE(j}Bi(gi - a 2 )||)) ' 



m,n=0 \ Z=0 

(Ag||g 2 ( a ; 1 -x 2 )||) m+w - ia (A/?) m+ " \ 

(m-/)!(n-Z)! j ' 1 J 

and 

VAR(V k )= [ Cov {V k { Xl ),V k (x 2 )) d Xl dx 2 . (8) 



2.2 Limit Laws 



We are now ready to obtain the limit laws by letting A — > oo and scaling the shapes by 5. Let C ((5, A) 
be the PB model C in which the shapes are scaled by 5, i.e., the shapes have the same distribution 
as SC. Let V k = Vk(X,S) be the resulting k- vacancy in R. 

Theorem 1. If 5 — > as A — > oo suc/t £/iat e) d A — > p /or < p < cc in the scaled coverage process 
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C(5, X), then 

a.s. 



3=0 J 



Theorem 2. Consider the scaled coverage process C(S,X). If S — > as A ^ oo sitc/i i/tat 5 rf A — ► p 
where < p < oo, then 



m) - wie , . ( 9 ) 



mi 

m=0 



EdVit - E(y fc ) H - /orl<p<oo, (10) 

AVAR(Vfc) -> a 2 , (11) 



where 




m,n=0" B °' 2t \ Z=0 

(12) 



2! (m-l)\(n-l)l 



mini 



Theorem 3. // 5 — > as A —> oo such that 5 d X — > p /or < p < oo in the scaled coverage process 
C(5,X), then 

VX(V k -E(V k ))^ N(0,a 2 ), 



in 



distribution where a 2 is as defined in ilty) . 



It follows from Theorem [T] that if 5 rf A — ► oo, then — ► almost surely. However, note that 
Theorems [1] and [2] do not guarantee complete coverage with high probability for large enough A. 
We now consider such a requirement, key to which is the inequality (|13p below. 

For the following two theorems, we assume that d = 2, the operational area R to be the unit square 
[0, l] 2 , and the sensing areas to be discs of radius r > satisfying irr 2 < 1. This last requirement is 
only to give a compact expression in the inequality below. Since our interest is in the asymptotic 
behavior of the coverage process as A — * oo and r = r\ —* 0, this is satisfied for all large enough A. 
For any A, r > define the event Z\[r) := {T4(A, r) > 0}. For radii r\ > which are decreasing in 
A to zero, we shall abbreviate Z\(r\) by Z\. We have the following inequality. 

1 < Pr (Z x (r)) < 2e~^ 2 (l + A V(l + £ , (13) 



1 + _ v v " _ \ v Xitr' ^ i 



7 



where 

e _ 4(fc + l)! 

e- A7rr2 A(Avrr 2 ) fe ' 

A formal proof of (|13p is given in the appendix. The above inequality is an extension of Theorem 3.11 
of [T] for k > 1. A similar inequality is derived in [TU] (proof of Theorem 1) under the condition 
that the operational area R is a square of side length I — > oo, with intensity A = A(^) — > oo (see ((TJ) 
and the sensing area of the sensors is one, i.e., 7rr 2 = 1. The following result follows immediately 
from (fT3l. 



Theorem 4. 5uppose t/iat d = 2, and let r\ = log X+k lo ^ l ° g A+c(A) . As A -> oo, Pr(Z A ) -► if 
c(A) — ► oo. Further, if c(A) -»c£S, then Pr (Z A ) > ci for some constant c\ G (0, 1). 

Theorem H] gives the critical radius required for complete A>coverage with probability approaching 
1 in two dimensions. We now show that by taking the radius to be a bit larger than that obtained 
from Theorem [31 we can get a stronger and more stable complete coverage regime. The following 
discussion will become easier if we assume A = n where n is an integer. To make the above notion 
of strong complete coverage more precise, define the critical radius for complete coverage as 

r* := inf{r n > : V k , s (n, r n ) = 0}, (14) 

where Vk )S is the A;-vacancy in the unit square. 

Theorem 5. Let d = 2, and let Vfc ;S be the vacancy in the unit square. Let r* be as defined above. 
Then, almost surely, 

lim = 1. 15 

n^oo log n + k log log n 



Remark: Let < e < 1. The above result implies that by taking the radius 

^ = (l + £ ) l0gW + fcl0gl0gn , (16) 

the unit square will be almost surely, completely fc-covered for all n large enough. Thus, if n is 
large, by taking the above r n , which is eventually larger than the one given in Theorem HI we can 
ensure a complete fc-coverage regime that will not see vacancies even if the number of sensors is 
increased (with corresponding decrease in r n ). Further, the above result gives a strong threshold 
in the sense that if 

2 Jog n + A; log log n 

r n = (!-e) , (17) 

7rn 

then the unit square will not be completely A;-covered for all large enough n, almost surely. 
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2.3 Proofs 



Proof of Theorem [JJ To obtain the required result, first note that for the A;- vacancy in a unit 
cube D, Vfc(-D), the expectation is given by E(Vk(D)) = e^^Yl'jZo ft • Now observe that the 
two scaling regimes — (1) R is fixed with X5 d — > p and (2) 5 = 1, Ri = IR, and I — > oo with P being 
a Poisson point process of intensity p are equivalent. See Section 3.4 of pQ for more discussion on 
this. Theorem Q] follows using the same steps as in the proof of Theorem 3.6 in [I]. □ 

Proof of Theorem [2} Since the sequence {Vk(5, \)}s>o,x>o is bounded, ([9]) and (fTU|) follow from 
Theorem [TJ 



(fTTj) follows from (|7|), (J8J) if A J fixfi Cov (xm( x i)> Xn{x2)) dx\ dx2 converges to the product of the 
integral on the r.h.s of (|12p and (p||i?||e _2p ^) . This is shown below. 

Let B[ (y) = (y + 6C) n 5C and = (y + 5C) n 5C C . Making first the change of variable x = x\ 

and y = xi — X2 in ([7]) and then y to <5y, we get 



A / / Cov(x m (x),x„(a;2))(i^2rfa; 

J_R Jx-R 



= x[ [ e - 2Udp J2 ^ii^wii)^l^?iMJ!))! : 



m,n— \ /— 



iRJx-R 

(AE(||Sf(t/)||)) m+ "- 2 ' (X(3) m+n 



(m — l)\(n — 1)1 mini 



dy dx 



\ I I e -2xs d e v f e ^ d K\\^(s- l v)\\)y^ (AgE(||gi(g ^)||)) 

J R Jx-R m ,n=0 \ 1=0 ^ 

(XS d E(\\B 2 (5- 1 y)\\)) m+n - 21 (XS d 0) m+n \ 

dy dx 



(m — l)l(n — 1)1 mini 



= X6 d e- 2Sd ^ 



R JS- 1 {x-R) 



j« ml B livm ^ (^XE(\\ Bl ( y )\\)) 1 ^ (^AE(||B 2 ( y )||)) ro+ " V 



/! 



x 



(m-Z)!(n-0! 



nun 



(18) 
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Figure 1: Figure shows the sets A±, A 2 and ^3. The squares are shaded light gray. The squares 
are of length rr5 and the 'spacings' between the sqaures are of width 2rd. The light gray shaded 
region within R is A\, the darker region within R is A2 and the unshaded region within R is ^3 
(Adapted from pQ). 

where 

TB{ ) ' J 5 -H x -R)nB , 2T V h U {m -l)\(n-l)\ 

(5 d X(3) m+n \ , s 

- 7, dy. 19 
mini I 

fs(x) as well as the integrand in (|19p are uniformly bounded, since for any e > and all A > 
sufficiently large, we have <5 -1 (x — R) D i?o,2r C Bq^t, |<5 d A — p| < e and -Bj C -Bo,r- Therefore, from 
(|18|) . (|19|) and the dominated convergence theorem, we have 



E>Q,2t 



1=0 



(20) 



and 



A / Cov (xm^Oxnfe)) dx 1 dx 2 ^pe 2/9/3 / c dx = pe 2p/3 c ||-R||, 
as A — > 00. This proves (fTTj) . □ 

Proof of Theorem [3} Let r be a large positive constant. Divide all of ^k d into a regular lattice 
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of ci-dimensional cubes of side length (rrS), with each cube separated from its adjacent cubes by 
a (i-dimensional 'spacing' of width 2t5. Let A\ denote the union of those cubes which are wholly 
within R, A 2 the union of the spacings that are wholly contained in R and A3 the intersection of 
R with all those cubes and spacings that are only partially within R. Fig. [1] illustrates these sets 
for d = 2. Since Ai, A2 and ^3 form a partition of R, the vacancy within R may be written as, 

where Vu is the k- vacancy within the region Ai. Under the assumptions of the theorem, as 5 — > 0, 
the cubes get finer. Further the number of spacings is less than \\R\\ /(rrS) . Since the volume of 
the 'spacings' is (rrS)^ 1 x (2t6), we have 

||A 3 ||->0, and ||A 2 || < 2||12||r~ 1 . (21) 

From ©, we get 

fc-i 



VAR(V fc (i) ) = Y, f dx i f Cov( Xm (xi)Xn(x2))dx 2 

m,n=0 Ai Ai 
fc— 1 

= 5 d J2 I e-^fs^dxtKMd^WAil 

n J Ai 



(22) 



m,n=0 

where the last inequality follows from (|20p . From the inequality above and (|2ip . we get 



From d2U) and ([22]>. we also get 



lim AVAR(V, (3) ) = 0. (23) 

A— >oo V / 



lim lim sup AVAR (F fc (2) ) = 0. (24) 



A— >oo 



Thus, for any e > 0, we can choose r large enough so that AVAR(V fc (2) ) < e, for all A sufficiently 
large. Hence in order to obtain the central limit theorem that we are seeking, we need to concentrate 
on V^ 1 ^ and obtain a central limit theorem for it. Since AVAR(T4) — ► a 2 from (jlip . we need to 
show the following. 



VP-EVP 



^ k ' d ^iV(0,l), lim limsup|A(VARfy fc (1) ) -VAR(Vfc))| =0. (25) 



VAR(F fc 



(1) \ \ 1/2 r ^°° X^oo 



Let n = n(\) denote the number of cubes of side-length rr5 in A\, and let Di denote the z-th such 
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cube, 1 < i < n. Denoting the /c-vacancy in Di by Ui, we have vjp = ^ Ui. Under the scaling 
regime, each shape is contained in B$ T $ and the spacing between the cubes is 2t5, and so no shape 
can intersect more than one cube. Hence, given A, the U-s are independently distributed and we 
have 

VARjV^) = VAR(£/i) = nVAR(E/i) . 

i 

Let D be a d-dimensional cube of side rr with the same orientation as D\. For any two real 
sequences a n , b n a n ~ b n implies that a n /b n — > 1 as n — > oo. From (J7J), (jSJ) we can write 

VAR(V fc (1) ) = ne~ 25dx/3 V f f ( e m('i-^+SS)nSS\\) y (*E(||(ai ~ ^ + gg) D <5S||))' 

x (AE(||(x 1 -x 2 + ^)n(^) c ||)r + "- 2 ' _ (^A/?) m+ " \ 
(m — l)\(n — 1)1 mini J 

k—l r r / mAn fnCf] \ D I „fzl\\T3tt\\m+n-2l 



f»e-»* £ / / f e*^ £ (P E (ll^ilD) x (PE(||B 



m+n \ 

, , - ) dxidx 2 . (26) 
m!n! / 

The rest of the proof follows as in ([1] pp. 157-158). □ 

Proof of Theorem [5j Recall that Z n (r n ) = {Vk(n,r n ) > 0}. Suppose we show that for the choice 
of r n as in (|16p . we have Pr {Z n {r n ) infinitely often ) = 0, then it follows that 

hmsup < 1. (27) 

n~>oo log n + k log log n 

On the other hand, for the choice of r n as in (|17|) . if we show that Pr (Z^(r n ) infinitely often ) = 0, 
then we can conclude that 

ton inf T^fSL > i. (28) 

n— »oc log n + k log log n 

First we show (|2T[) . Take subsequence rij = j a , where a > will be chosen appropriately later. 
Define the events Ej = U™L%Z n (r n ) C Z n .(r n . +1 ). Then, Pr (Ej) < Pr {Z n .(r n . +1 )). 

For r n as in (|16p . observe that njr^. +i — ► oo. Hence, from (|13fl . we get, for large enough j, 

Pr{Z nj (r nj+1 )) < Cike- n ^ rl ^ nj (n^rl ]+1 )\ 
where Ci is some constant. Choose n > (1 + e)^ 1 . Since nj/nj+i f 1, as j — > oo, we have for all 



11 (m - l)l(n - l)\ 
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large enough j, rij/rij+i > rj. Hence, r/(l + e) logn J+ i < njirr^. +i < 41ogn J+ i, so that 

Pr [Z nj {r nj+1 )) < C2e-^ 1+e)losn ^(logn j+1 ) k nj < C 3 (log j) k (j + 

where C3 are constants that depend on a, k, e. Thus for a large enough we get Yl'jLi P r (-%) < 00 • 
By the Borel-Cantelli Lemma, Ej happens only finitely often almost surely. By the definition of 
Ej, this implies that Pr(Z n (r n ) infinitely often ) = 0. 

Similarly using the lower bound for Pr(Z n ), we can prove (|28p . As above, let rij = j a , j > 1 and 
a > to be chosen later. Let r n be as in (fT7|) . Using the lower bound in (fl~3j) . we get 

Pr (U^,^(r„)) < Pr (z c n . +1 (r n .)) < < 6 3 , 

where 

4(fc + l)! 



e Jn j+ i(n j+ i7rr^) fc 

Proceeding as in the proof of the upper bound, we can show that the above probability is summable. 
This completes the proof of the theorem. □ 



3 Path Coverage in the Markov-Poisson-Boolean Model 
3.1 Preliminaries 

Let C t = {Y{(Xi +Ci), i > 1} be the MPB model as defined in SectionO Recall that {Y^t > 0} 
are i.i.d. copies of {Y l ,t > 0} which is a {0, l}-valued continuous time stationary Markov process. 
The sets C{ are independent and distributed as C satisfying C C Bq :T for some fixed r > 0. Let /j,q 
be the transition rate from the 0-state to the 1-state and fj,\ the transition rate from 1-state to the 0- 
state. This of course means that in each visit, Y is in the and 1 states for exponentially distributed 
times with parameters (io and \x\ respectively. Then, (see [6], Chapter 6), the stationary probability 
of the sensor being in state j is given by pj := Pr (Y t = j) = ; j S {0, 1}. The transition 

probabilities between the states are defined by pt(j, k) := Pr (Y s+t = k\Y s = j\ for j, k S {0, 1}. It 
can be shown that 

Pt(j,j) = (1 -Pj)e" 7 * +pj, 
where 7 = ^ + W and p t (j, k) = 1 - p t (j,j) for j, k £ {0, 1} and j / k. 

Now, consider a target moving on a straight line path in 3? rf of length L units with a velocity of c 
units per second. As in the previous section, we will let X{ = (Xf,Xf, . . . ,Xf) be a point in 3f? d . 
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Without loss of generality, we can consider the target to be moving along a coordinate axis, say X , 
for T := L/c units of time. Define s := (cs, 0, . . . , 0) G K 2 , where s G s is the position of the 
target on the X-axis at time s. Let be the indicator variable for the target being covered by 
one or more of the sets Cj at time s, and = V\ t be the random variable denoting the duration 
for which the target is not covered in [0, T], i.e., 

V T ■= / T (l - CM) (29) 
•/ o 

The effect of the transitions of makes the study of Vr interesting. We first calculate the 
expectation and variance of Vr- Since the Y? are stationary, the expectation of Vr is straightforward 
and is given by 

E(V T )= [ T Pr(s(jtC s )ds = Te- Xpif3 , (30) 
Jo 

where (3 = E(||C||) . The second moment of Vr can be obtained from 

E[v£) = [ [ Prfei £C sl ,s 2 <£C S2 ) d Sl ds 2 
Jo Jo 

= 2 / Pr(s! <£C sl ,s 2 (£C S2 ) dsx ds 2 . 

Jo J si 



We now evaluate this integral. Define ^i(si,s 2 ) := (s t + C) D (s 2 + C), B 2 (s 1 , s 2 ) := (sj + C) n 
(s 2 + C) c and -Bs(s 1 , s 2 ) := (s 2 + C) n (s 1 + C) c . From this definition, B\, B 2 and S3 are disjoint 
and we can see that H-B2II and 1 1 .S3 1 1 have the same distribution and ||-Bi|| = || — s 2 + C) D C||. 
We will suppress the obvious arguments of the sets Bi. As in the arguments leading to ([6]), we 
consider a set A' = B Sli2r U B s ^ )2t and assume there are iV sensors within A'. Then the integrand 
in the above equation is 

Pr (0 on sensors covering s x at si; on sensors covering s 2 at s 2 ) , 

which equals Pr (Mi = 0, M 2 = 0, M3 = 0, M4 = N), where Mi = no: of sensors active at s± or s 2 
among those in A' which can cover s x and s 2 , M 2 = no: of sensors active at s\ among those in A 1 
which can cover s x and not s 2 , M3 = no: of sensors active at s 2 among those in A' which can cover 
s 2 and not s l3 M4 = remaining sensors in A'. Now note that 

(Mi, M 2 , M3, M4) ~ Multinomial(A^, 01, 02, 03, 04 = 1 — a% — a 2 — 03), 
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where 



a x = (l-p p aa - Sl (0,0))E([|5i||)/[|A , [| J 
a 2 = piEGIBalD/M'H, 
a 3 = piEdl^lD/KH. 

a\ is the probability that a sensor can cover both of s x and s 2 an d is active at either si or s 2 - It 
is the product of probabilities of a sensor lying in B\ and making the desired transition of states. 
The first probability is E(||i?x||) /||A'||. The second probability is obtained as follows: 

Pr (a given sensor is on at s\ or s 2 ) = Pr (the sensor is on at s\) 

+Pr (the sensor is off at s\ and on at s 2 ) 
= Pi +P0Ps 2 -si(0, 1) = 1 -poPfl 2 -si( > )- 

The remaining ctj's are even simpler to obtain by similar calculations. Since N follows a Poisson 
distribution with mean A||^4'||, the integrand in the above expression for E(V^) is given by 



00 
£=0 



l ('MI^'ll) P -A||A'|| _ ^-ACai+as+aajm'l] 



= 2 [ T /' T e -2piA/3 e -AE(||B 1 ||)(l-2p 1 -poP S2 - sl (0,0)) ds2 rfsi _ 



Hence, 



VAR(y r ) = E(y T 2 ) - (E(V T )f 

= 2e- 2 ^^ ^ ^ ^ e -A^||flx||)(i-Pop^-. 1 (o,o)-2 P i) _ ^ ds2 dsv (31) 

Observe from (I30p that the expected coverage depends neither on the shapes of the sensing regions 
nor on the transition probabilities of the on-off process. Their effects are picked up only in the 
variance of the coverage as can be seen in (|3ip . 
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3.2 Limit Laws 



We can now delineate the scaling laws. Let Cf(S, A) be the process Ct in which the shapes C are 
scaled by 5. In addition, the parameters of the Markov chain governing the on-off process will also be 
made to depend on 5. Without loss of generality, let fJ-o(,$), fJ-l(8) be such that Ho(5) / Hi(S) = Ho/pi 
for all 5. It can be easily seen that the proofs remain valid when po(5) / pi(5) — > po/pi as 5 — > 0. 
Hence the scaled stationary probability is Pi(5) = pi and the scaled transition probabilities are 
pf(0, 0) = pie~^ s ' s + po, where j(5) = /j,q(5) + fJ>i(5). Observe that the case //o(5)///i(<5) —> 0(oo) 
imply that po = 1(0) and they correspond to sensors being always off (on). 

We will first state the four main theorems for the MPB model and provide their proofs after a brief 
discussion. 

Theorem 6. Consider the scaled coverage process Ct(5, A). Let 5 — > as A — > oo such that 5 d X — ► p 
where < p < oo and 5~f(5) — > ao> where < oq < oo. Then, 

E(V T ) -> Te- pi ^, (32) 
E(| V7 1 - E(Vr) | p ) -> /or 1 < p < oo. (33) 

^VAR^t^)) -»■ (7?(ao), (34) 

where 

o-f(ao) = 2Te- 2pi ^ ^~ (>^ll(s+C)nq|)( w+ poe-«o*) _ ^ ds (35) 
Theorem 7. Let i/te assumptions be as in Theorem® Let < ao < oo. Then, 

r*(Vr(6) - E(V T (5))) - JV(0,<^), (36) 

in distribution, as 5 — » 0. 

The following two theorems describe the asymptotic behavior of the coverage process Vr as the size 
of the operational area becomes very large. 

Theorem 8. Consider the coverage process Ct = Ct(l, A). As T — » oo, 

V T e~P lX(3 

— > a.s. 

T c 

Theorem 9. For the coverage C t = C t (l, A), as T — > oo, 

T- l '\V T - E(V T )) -> iV(0, ffl), (37) 
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in distribution, where 



■It 



a 2 = 2e~ 2piXf3 / C ( e Pi AB ;ilfe+ c, ) nC 'll)(Pi+Poe- (A10+Ml)3 ) _ i ) d s 



Remark: The scaling of the shapes by 5 can also be viewed as follows: If the velocity of the target 
increases by 1/6, then the time it shall spend in the region of a sensor shrinks by a factor 5 but it 
shall travel a larger distance and hence be seen by more sensors. In light of this remark, one can 
view the scaling results in the unreliable sensor networks case as tracking of a high- velocity target 
in a highly fluctuating sensor network. Observe that 

o?(0) = 2Te~ 2 ^ f 2 j' C (>pB[ll(i+C)nC||) _ ^ ds _ 

Since the target is moving fast, the contribution to the variance is coming from locations that are 
close to each other. When ao is small, the target moves much faster as compared to the rates of 
transition of the sensors between the on and the off states. Thus given that there was a sensor 
on at a given location, the chance of it being on in the near future will be close to 1. Similarly if 
ao = oo, then 

<x?(oo) = 2Te~ 2 ^ £ T/C ^ipKWis+QnCW) _ ^ dfj 

Hence when ao is large, the transitions happen at a much faster rate compared to the speed of the 
target. So given that there was a sensor on at a location, the chance of the target seeing a sensor 
on in the near future will be roughly the stationary probability of a sensor being on, which is p\. 



3.3 Proofs 

Proof of Theorem [6} ()32p is straightforward. For (|33p . as in the proof of (jlpp . it suffices to show 
that the variance converges to 0. Putting s = S2 — si in (|3"T]) and the expression for p S2 _ Sl (0,0) 
and simplifying, we get 

VAR(V T (5)) = 2e~ 2Sdx ^ £ £~ S1 ^pA^^nc^+poe-^) _ ^ ds d ^ 
Replacing s by s5 we get 

T S — ^ (T s ^ 

VAR(V T (5))=25e- 2Sdx ^ J ^ P ^||(, + C)nc||)( Pl+P o e -^» s ) _ ^ ds dsv (38) 

The range of integration in the inner integral diverges to oo as 5 — > 0, and the limit of the integrand 
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is as described in (|35p (using the fact that C C Bq <t ). We can justify the convergence to limits since 
the inner integral and the integrand are uniformly bounded and hence the dominated convergence 
theorem applies. Indeed, for any e > fixed and sufficiently large A we have X5 d < p + e and 



□ 



Proof of Theorem We shall first prove (|36p . Proof follows the same idea used for proving the 
central limit theorem in case of /c-coverage. 

Choose a large constant r. Divide [0, T] into alternating intervals of lengths re)r (type 1) and 
25t/c (type 2), where c is the velocity of the target. Truncate the interval containing T at T. 
Denote the union of type 1 intervals by A\ and the union of type 2 by A<i- Let the vacancies 
arising in A 1 be denoted by and that in A 2 by V^ 2 \ Now, V T = + W 2) . Note that 
1 1 ^2 1 1 < (28- x -?-) = ^— >0asr— >oo. In ()38|) . the inner integral converges and hence it is 
bounded. Using this fact when the range of integration is Aj, we get 

VAR(V W ) < C5\\Ai\\. 

It follows that 

var(t/( 2 )) , 

lim limsup = 0. (39) 

Thus we need to show the following to prove the central limit theorem. 

i =-4 JV(0, 1), lim limsup |A(VAR(F (1) -VAR(y T ))| -» 0. (40) 

(VAR(y( 1 ))) 1/2 



r^oo 



A— >oo 



Let n = n{8) denote the number of intervals in A\. Let Di denote the i-th interval of type 1. Let 
Ui be the vacancy in Dj. We can see that the Ui are i.i.d. random variables except for last one 
corresponding to the truncated interval and = Y12=l ^i- As in ([38|) . we obtain, 

VAR(V (1) ) ~ nVAR([/i) 

= n25e- 2Sdx ^ j* ^ d ^iB(ll(,+C)nc||)( Pl+P oe-^w=) _ ^ ds ^ 

~ n5 2 2e~ 2 ^ J J ^ e PPiB[||(i+a)na||)( w + OT e-o^ _ ^ dg ^ 
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The last relation is obtained by replacing s\ by 5s\. Since E(|f7j — E([Zj) | 3 ) < ||D^||VAR(C/j) = 
r5TVAR(Ui) and n = 0(^ 1 ), as 5 -> 0, 

Z,W,-WJl? rjr = , ^ 

(EiVAR(t/,)) 3/2 (E^ARfC/,)) 1 ' 2 



Hence the first part of (|40j) follows from Lyapunov's central limit theorem and the second part of 
the proof follows as in p] (pp. 157-158). □ 

Proof of Theorem [8} We divide the interval [0, T] into intervals of length 2r/c where the number 
of odd numbered intervals is n and the number of even numbered intervals is m. n and m are such 
that n + m = [cT/2t\ and < n — m < 1. Let U} be the vacancy of the i-th odd-numbered interval 
and Uf the vacancy of the i-th even-numbered interval. Hence, 

cT 



V T = I l[s(£C s ]ds 
Jo 

n »2(2i+l)r/c m /•4(i+l)r/c /-T 

= J2 l[s^C s ]ds + y2 l[siC s ]ds+ l\s£C 

,_n J4jr/c ■ n ,/2(2i+l)77c J(n+m)2r/c 



i=0 

i=0 i=0 i(n+m)2r/c 



As T — > oo, 



1 /" T 2r 
- / l[s ^C s ]<fe < — ^0. 

1 J(n+m)2T/c 01 



' {n+m)2T/c 

{J7j } and {t/f} are sequences of i.i.d. random variables while \J\ and Uf might be dependent for 
some Vs. Hence, by the strong law of large numbers, 



n 

-1 ' 

n 

i=0 



1 ^U^E(ul)=2re-^/c a.s., 



and 



Since 



m 

—i 

m 

i=o 



^U^E(u!)=2re-^/c a.s.. 



T r n m 

Vt ( n 

' n 



W/ \^ J J(n+m)2r/c 



T \7\ 

the strong law follows by noting that n/T,m/T — ► l/4r as T — > oo. □ 

Proof of Theorem [9l The proof follows using a similar spatial blocking technique as in Theorem 
[7] provided we show T _1 VAR(Vj') — > er| as T ^ oo. First make the change of the variable s = S2 — s± 
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and then u = s\/T in (|3ip . to obtain 

VAR(Vr) = 2Te- 2p ^ C /™ (^mts.+^nc^-^e-^^) _ 1 j dg du _ 



o Jo 



The result now follows from the monotone convergence theorem and using the fact that C C -Bo,t- 
The finiteness of the limiting variance follows from uniform boundedness of inner integral in the 
above equation. Indeed for sufficiently large T, we have 

T(l—u) r-T 

' , P iAB;||(s+C)nC||)(pi+poe-^o+Mi) S ) _ ^ dg < j ^gpiA^llfe+CTJnOIDCw-Hwe-Cw+M)*) _ {\ dg 



.piAEdKs+^nCIDCpi+poe-^o+fi)^ 







< ^ fallen) + ^ <QO _ 



c 



□ 



4 Appendix 

Proof of Inequality I13t Let R be the unit square [0, l] 2 . Consider the PB model defined 
in Subsection 11.21 with intensity A and the sensing regions to be discs of radius r > 0, satisfying 
7rr 2 < 1. Let (see Q) be the k— vacancy in the region R. First we shall derive the upper bound. 
We can write 

Pr(V k >0) =Pi+P2+P3, (42) 

where 

pi = Pr(no disk is centered within R) = e _A < e _A7rr , since irr 2 < 1, (43) 
P2 = Pr (at least one disk is centered within R but none of the disks intersect any other disk and 
none of the disks intersect the boundary of R) 
< Pr (at least one disk is centered within R) x Pr (a given disk intersects no other disk) 
= (1 - e~ x ) x e - Xwi2r)2 < e" W , and (44) 
P3 = Pr (R is not completely A:-covered, at least one disk is centered within R and 

at least two disks intersect each other or at least one disk intersects the boundary of R) . 

We shall now obtain an upper bound for p% using the same technique as in [TO]. Define a crossing 
to be either the point of intersection of two disks or the point of intersection of a disk with the 
boundary of the region R. A crossing is said to be /c-covered if it is an interior point of at least 
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k disks. By Theorem 4 of [9j, R is completely A;-covered iff there exist crossing points and every 
crossing point is /c-covered. 

Let N and M be the number of crossings and the number of crossings that are not /c-covered 
respectively. Then, by definition of crossings and its relation to complete /c-coverage, we have 

P3 < Pr(M>2)<^ 

< — E(JV) Pr (a crossing is not fc-covered) 

k ~ l /a r 2\f 
1 1 



i(E(^) + E(N 2 )) (e-^ 2 £ , (45) 



\ i=0 

where N\ is the number of crossings created by intersection of two disks and N2 is the number of 
crossing created by intersection of the disks with the boundary of R. By Palm theory for Poisson 
processes (see [8]), we have the following 

E(iVi) = E(No. of pts in R) E(No. of crossings created by a given node with other nodes) = 4A 2 7rr 2 . 

(46) 

In the above formula, the two crossings created by intersection of two disks are assigned one each to 
the two nodes. For estimating E(^) note that only nodes within a distance of r from the boundary 
can create crossings with the boundary. The expected number of such nodes is at most 4Ar and 
each node can create atmost 2 crossings. Again by Palm theory, we have 

E(iV 2 ) < 8Ar. (47) 

Substituting from ([56]) . (|T7]) in (jl5) . we get 



k-l 



2 \„-A7rr 2 V- ( X ^ r 



P3<2A 2 vrr 2 (l + -^)e- W ^ 



Xirr ^— ' i\ 

i=0 

The upper bound in now follows from the above inequality and (|42|) - (|44|) . The lower bound is 
based on the following inequality which follows from (3.7) of pQ. 



(E(^)) 2 

ri s U) / - 

By definition 



Pr(y fc >0)> VSf- (48) 



E(^ 



E(y fc 2 ) =h + h, 
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where 



h = I I E(V k (xi)V k (x 2 )) dxt dx 2 < (E(V k )) 2 , 

lR 2 n{\x 1 -x 2 \>2r} 



h = E(V k (xi)V k (x 2 )) dxi dx 2 . 

J J R 2 n{\x 1 -x 2 \<2r} 

The bound for I\ above is obtained using the independence of the two terms in the integrand. Let 
B Xi ,r be disks of radius r centered at Xi, i = 1,2 respectively. By standard calculations (see proof 
of Theorem 3.11, pQ), it can be shown that ||-Bx 2 ,r \ Aci^H > irr\xi — x 2 \/2. 

E{V k {xi)V k {x 2 )) dxi dx 2 

R 2 n{|xi-x 2 |<2r} 

/ / Pr(both B xi ^ r and B X2 ^ r contain less than k nodes) dx\ dx 2 

J J R 2 n{\x 1 -x 2 \<2r} 



< / / Pr ( B Xltr contains less than k nodes, B X2 ^ r \ B Xl<r contains less than k nodes) dx\ dx 2 

J J R 2 C\{\x 1 -x 2 \<2r} 

= P r {B Xl ,r contains less than k nodes) Pr (B X2 , r \B Xl , r contains less than k nodes) dx\dx 2 

J J R 2 C\{\xi-x 2 \<2r\ 



< 



R 2 n{\x!-x 2 \<2r} 

E(Vfc(xi)) (e-^-x^y i^-x^Y . 

R?n{\ Xl -x 2 \<2r} V ll J 



E(y k (xi)) dxx / e-wzy \«>"h*) 27rxdx 

R Jo rr? i! 



nr 2 



x-ktx/2 (A7rrx/2) 

i- 

Jo f=i i! 



(Vfc 



i=0 

fc-1 



JU 2=0 



where the equality in the penultimate line is obtained by making the change of variable u = irrx/2. 
The last equality follows from the usual representation of the Gamma function (see (26), [10J). 
Substituting the above bounds for I\ and I 2 in ([^8]) . we get 



1 + x 2 -Kr 2 By k ) 

Setting (3 = irr 2 in ([5|) and ignoring all but the last term, we get 

"■'(Avrr 2 



(49) 



-Xirr 1 2\k-l 
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The lower bound in (fT3|) now follows from (f4T)|) . ([50]) . □ 
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